Domain walls in three-field models 
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We investigate the presence of domain walls in models described by three real scalar fields. We 
search for stable defect structures which minimize the energy of the static field configurations. We 
work out explict orbits in field space and find several analytical solutions in all BPS sectors, some 
of them presenting internal structure such as the appearence of defects inside defects. We point out 
explicit applications in high energy physics, and in other branches of nonlinear science. 
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This paper deals with domain walls presenting inter- 
nal structure that may engender nontrivial phenomena. 
Specifically we propose new systems of three coupled real 
scalar fields whose dynamics are controlled by superpo- 
tentials in the bosonic sector of supersymmetric theories. 
The systems allow for a multiplicity of solitonic solutions 
including the presence of defects inside defects, introduc- 
ing new results that lead to a variety of applications in 
diverse fields of physics, chemistry and biology. 

Domain walls may involve energy scales as different as 
those found in fields of magnetic systems nonequilib- 
rium thermodynamics |^ and cosmology Usually, do- 
main walls are immersions in (3, 1) dimensions of linear 
kink-like defects living in (1, 1) dimensions. These ex- 
tended spatial structures are the basic building blocks to 
describe the coexistence of different domains or phases, 
and they are usually modeled by a set of scalar fields 
playing the role of order parameters satisfying amplitude 
equations. They also appear in models that include su- 
persymmetry [|| ^, ^, 0| and supergravity ||, ||] . In par- 
ticular, a great deal of attention has been given to the 
restrictions imposed by supersymmetry in the classifica- 
tion and solution of the solitonic field equations - in a 
supersymmetric theory one can categorize the topologi- 
cal defects as BPS and non-BPS states, according to the 
work of Bogomol'nyi and of Prasad and Sommerfield pO[ . 

In models described by a single real scalar field we 
may distinguish at least two classes of systems, one that 
supports a single type of defect, and the other, that sup- 
ports two or more different defects. We exemplify this 
by recalling the (/)^ model, which supports a single type 
of defect, the tanh-like kink, and the double sine-Gordon 
model, which may support two different defects, the large 
and the small kinks ||ll| . On the other hand, models con- 
taining two or more real scalar fields give rise to at least 
two other classes of systems - those_that support defects 




and, those that 
. In the latter con- 



that engender internal structure 
support junctions of defects 

text, in Ref. jwj, one has investigated the possibility of a 
regular hexagonal network of defects to spring in a model 
described by two real scalar fields. In Ref. |15 the idea 



recently presented in 
that entraps a fulerene 



of an hexagonal network of defect to be nested inside a 



topological defect has been investigated in a model de- 
scribed by three real scalar fields. That investigation 
has shown that when the host domain wall is driven to 
relax to cylindrical or spherical shape, the nested hexago- 
nal network could give rise to structures resembling nan- 
otubes or fulerenes, respectively. Similar ideas have been 
1^, in the case of a soliton star 
ike network of domain walls, and 
in |17| , where the Skirme model is used to show the pres- 
ence of fulerenelike structures in light nuclei, and also in 
p8[ , in the case of solutions of the Einstein equations 
describing black holes pierced by cosmic strings, forming 
structures that very much resemble the Platonic solids. 

In Ref s. [pi], one has given up supersymmetry, to cir- 
cumvent issues concerning stability of the nested network 
of domain walls. Here, however, we return to the basic 
problem, which concerns investigating models described 
by three real scalar fields in the bosonic portion of super- 
symmetric theories. Models with three real scalar fields 
have been studied before, for instance in an 0(3) model 
which admits unstable embedded walls jl^, and also in 
pO| , which investigates stability of kinks in a deformed 
0(3) linear sigma model. Other motivations include the 
possibility of building spatial junctions of defects, instead 
of the planar junctions considered in ^,0, and the issue 
of finding a host domain wall that entraps a network of 
defects within the framework of supersymmetry. To con- 
struct spatial junctions of defects one first recalls the sim- 
pler issue, of contructing planar junctions. As we know, 
we can tile the plane with regular polygons in three dis- 
tinct ways: with triangles, squares or hexagons. These 
cases require junctions with six, four and three legs, re- 
spectively. Thus, the two-field model should engender 
discrete Zq, Z4, or Z3 symmetry. Interestingly, however, 
we notice that the triangular lattice is dual to the hexag- 
onal lattice, and vice versa, and that the squared lattice 
is self-dual. This means that there is a dual relation be- 
tween the (triangular, square or hexagonal) disposition 
of minima in the two-field model, and the (hexagonal, 
square or triangular) lattice made with the junction that 
spring with the very same minima. We take this as a 
guide to construct spatial junctions of defects in three- 
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field models, to fill space with regular polyhedra. Before 
investigating such issue, however, let us first deal with 
the more general problem of finding topological defects 
in three-field models engendering discrete symmetry. 

In general we are interested in potentials that can be 
obtained in terms of superpotentials, in the form 



2 2 2 



(1) 



where W^^ = dW/d(f>i, 4>i — (t>,X:P- The energy den- 
sity for static solutions that depend on the single spatial 
coordinate x can be written as 

i ^ ^ 

It can be minimized to give Eah — with 

AW'ab = Wa-W^. (3) 

where the indices identify vacuum states, such that 
= W^(0a,Xa,Pa) and SO forth. This energy is the 
energy of Bogomol'nyi-Prasad-Sommerfield (BPS) solu- 
tions The BPS states obey the first order ordi- 
nary differential equations 

^^-'^ w 

We now introduce the model, given by the superpoten- 
tial depending on two real parameters, r and s. 

Here we are using dimensionless space-time coordinates 
and fields. The corresponding potential is 

V = i(l-^2)2^2rW + i^^(/ + pY + 



2 „2 



(6) 



It may support two, four, and six minima, depending on 
the values of the two real parameters r ^ and s ^ 0. 
In the case of six minima, for r > and — 1 < s < 1 we 
get 



VU2 - (±1,0,0) 



V3,4 = I 0,±\/-,0 



V5,6 



(7) 

(8) 

(9) 



The model engenders the Z2 x Z2 symmetry, correspond- 
ing to reflections in the x and p axes. The presence of the 
s-dependent term in the superpotential forbids reflection 
symmetry in the (f> axis. However, we notice that the op- 
eration (f) —f —(j) leads the system to a partner system, in 
which one changes s — > — s. 



Before studying the topological solutions of this sys- 
tem, it is interesting to analyse the projections over di- 
verse axis and planes. The projection of the potential 
over the axis gives 



1 
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F(0,O,O) = -(l-0^) 



(10) 



which is known to produce kink-like domain walls that 
may host internal structure. In the (0, x) plane we have 

F(0,X,O) = i(l-02)2 + 2r20V 
We see that inside the domain wall 



v{o,x,o) = lr^(--r 

2 \r 



and that outside the wall 



V{±l,x,0)^2r\^ + -r\'' 



(12) 



(13) 



Thus, if the 4> field gives rise to a domain wall, then inside 
this wall the x field generates another defect. Notice that 
the squared masses of the x field inside and outside of the 
domain wall are given by m^(zn) = 4r and TO^(out) = 
4r^ , showing that the parameter r may be used to control 
the presence of elementary x mesons inside or outside the 
host domain wall, according to r being greater or smaller 
than unit, respectively. 

The projection of the potential in the (0, p) plane gives 

1 



^(-^,0,^) = -{l~^y + 2r'i^~sfp^ 
-r(l-02)p2 + irV 



(14) 



which takes the following forms, inside and outside the 
wall 

V{0, 0:P) = l+ r{2rs^ - l)p' + P^ (15) 



F(±l,0,p) = 2r2(lTs)V' 



^24 



(16) 



Clearly there is no spontaneous symmetry breaking 
(SSB) outside the host domain wall, but for 2rs^ < 1 
there is SSB inside the host domain wall and the sys- 
tem may generate domain walls with internal structure. 
The mass of the p meson inside the host wall, in this 
case, reads m^{in) = 4r(l — 2rs^). Outside the host wall 
there is an asymmetry in the -I- and — sides which is 
controlled by the vacua vi,2 — (±1,0,0) and the mass 
of the p meson is m'^{out,±) — 4r^(l ± s)^. The s pa- 
rameter induces the asymmetry for the p meson, since 
its mass depends on which side of the wall the meson is. 
Despite this asymmetry, the ratio rn^Jinj/rn^JoutjzL) = 



3 



(l/r)[(l-2rs2)/(lips)2j shows that for 1/r < l-2|s|+3s2 
the p mesons prefer to hve inside the host waU. We notice 
that the two conditions 2rs^ < 1 and 1/r < 1 — 2|s| + 3s^ 
define a large region in parameter space, and for s very 
smaU we have r > l+2|s|. The special case rs^ = 1 is also 
worth mentioning since it defines an interesting region in 
parameter space, in which there are explicit analytical so- 
lutions connecting the minima ^3.4 to v^, q by linear orbits 
(see below). In this case, as discussed above, the host wall 
cannot entrap topological defects but the masses of the p 
mesons can be explicitly computed to give m'^{in) — 
inside the host wall and mp(out,±) = 4(1 =F s)'^ / out- 
side. Thus, m'^p{in)/m'j,{out,±) = (l/2)[s/(l=Fs)]^, show- 
ing that, even without explicit SSB but with s sufficiently 
small, the p mesons prefer to live inside the host wall. 
The fact that p mesons prefer to live inside the wall is 
new, it does not appear in the other case, in the plane 
((/), x). More importantly, although s induces an asym- 
metry for the p field, we may control both r and s to 
entrap p mesons inside the host wall, in a way such that 
they do not perceive the asymmetry the model engenders. 

The model defined by Eq. (|^) has several topological 
sectors. All the sectors are BPS, except those connecting 
the minima ^3 to W4 and to vq , which are of the non 
BPS type. The energies of the BPS sectors are given as 
follows. In the vi^2 minima sector we have E]^pg = 4/3. 
In the sectors connecting vi^2 to 'i;3_4 we have Egpg = 
2/3. For the connections between wi,2 and v^^e 



E 



BPsi^) ~ 3 ^ ^ 



and for connections between U3 4 and W5 g 



E 



BPS 



is) 



(17) 



(18) 



We notice that these energies do not depend on r. They 
are plotted below, to show their behavior with the pa- 
rameter s. Eq. ( p^ ) shows that the limit s ^ turns 
to zero the energy in the sector connecting the minima 

This is so be- 



(0,±7l/r,0) and {s,0,±y/{l - s'^)/r). 
cause in the limit s — > these minima degenerate to a 
continuum set of minima, which may be connected with 
no energy cost, according to the Goldstone theorem. It 
is interesting to see that E^pg{s > 0) — Ebps{s > 0) 
and E]^pg{s < 0) — Ebps{s < 0) do not depend on 
s, and that there are two specific values of s [s = 
±s, with - 3s + 1 = 0, < s < 1] for which 
Espsi^s) — E'^pg{±s); for these two values the ener- 
gies of all the BPS sectors become equally spaced, or- 
dered as -Erp^^, E'rpc for s > or -Erpq for s < 
for 







^BPS^ 



and E 



BPS 



s > or 



'BPS 

E+ps for 



s < 0, 



with the values 4/3, 1, 2/3, 1/3, respectively. 

Let us now examine the presence of defects in the s ^ 
case. The BPS states obey the first order equations 
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FIG. 1: Energies of the BPS states given by Eqs. and 
( p^ as functions of the parameter s. 



dx 
dp 

dx 



-2r0x 
~2r(0 — s)p 



(20) 
(21) 



In the sector connecting the minima vi^2 there are topo- 
logical solutions that come from a straight line orbit with 
X = p = 0, which is described by 4>{x) = tanh(a;). This 
is the one-field kink solution for the system projected 
into the </> axis that generates the domain wall hosting 
the other fields, as we referred to above. For the system 
projected in the {4>, x) plane we get the superpotential 



VK((/>,x)=0-:70^-7-0x 



(22) 



This model was first examined in |2l|, and also in p2[ . 
In this case all the orbits in (0, x) plane can be obtained 
explicitly, as it was recently shown in Ref . . 

For the system projected in the (</>, p) plane the super- 
potential is 



w{4>,p) = 4>-\4>''-r{4>'s)p' 

It gives rise to the first order equations 



dx 
dp 

dx 



l-cj)^ - rp^ 



= —2r{(f) — s)p 



(23) 

(24) 
(25) 



To find solutions that connect the minima W1.2 to W5.6 
we follow the approach developed in Ref. . We found 
straight line orbits connecting these minima with 



;i 1 



(1- 



(26) 



and 



dx 



1 - - r(x' + p') 



(19) 



± 



1 1 



r 1 



(1 



(27) 
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For the first trajectories, Eq.(26), there are BPS states 
for \/r = 1 — s. They are 



{x) = i(l + s) + i(l-s)tanh(x) 
1 



(28) 



p[x) = ±-(l-s)Vr+s[l-tanh(x)] (29) 

They obey <j){x — > — oo) — s- s and p{x — > — oo) ^ (1 — 
s)Vl + s, and 0(a; — > cx)) ^ 1 and p{x — oo) — s- 0. 
Thus, they connect {s,0,±y/{l ~ s^){l - s)) to (1,0,0). 
For the second trajectories. Eg. (|27|) , there are BPS states 
for 1/r = 1 + s. They are 



p{x) 



-i(l-s) + i(l + s)tanh(a;) (30) 
1 



±-(l + s)Vl~^[l + tanh(a;)] (31) 



They connect (-1,0,0) to (s, 0, ± v^(l - s2)(l + s)). 

There are other orbits connecting the minima wi,2 and 
Vbfi- We found the following parabohc trajectories 



(32) 



giving rise to BPS states connecting the minima for r 
1/4. They have the form 



i(l + s) + i(l-s)tanh[x/L-(s)] (33) 



p(x) ^ ±[2(l-s2)]i/2{l-tanh[x/L-(s)]}i/2(34) 



and 

<P(x) = -i(l-,s) + i(l + s)tanh[a;/L+(s)] (35) 
p(x) = ±[2(l-s2)]i/2{l + tanh[x/L+(s)]}i/2(36) 

where L=^(s) = 2/(1 ±s) give the width of the topological 
defects. 

We also found three-field BPS states connecting the 
minima v^^^ to v^^q. This is even harder, because no 
field can vanish along any smooth orbit connecting 113.4 
to W5.6. We first notice that there are four different ways 
to connect (0,±yT7^,0) and (s,0,±V(l - s^)/r). We 
found straight line orbits described by 



' ± s-\/rx — 3 = 



and 



which imply that 



s^)(f>Tsp = 



(37) 



(38) 



(39) 



These orbits also project as straight line segments in each 
one of the three planes (</>, x), (0, p) and (x,p)- We use 
them to rewrite Eq. ([l^) in the form 



d0 
dx 



which is solved by 



(a;) ^ 2^'^ + tanh(a::/s)] 



The other solutions are 



Xix) = ±iy i[l - tanh(a:;/s)] 



(40) 



(41) 



(42) 



p{x) 



±iwi(l - s2)[i + tanh(x/s)] (43) 
2 V 7" 



They require that rs^ = 1. 

As we have shown, we have been able to find explicit 
BPS solutions in all the BPS sectors of a model described 
by three real scalar field with discrete symmetry. The 
presence of the parameter s restricts the symmetry of the 
model to Z2 x Z2, since s 7^ forbids reflection symmetry 
along the (j) axis. For this reason, now we study the 
behavior of the system in the limit s ^ 0. In this case 
the superpotential (||) becomes 



W = 



1 



^'-r^ix^+p') 



(44) 



The corresponding potential supports the two minima 
(±1,0,0), and also (0, ^/iJrcosO, y^l/r sin0), which de- 
fines a circle in the (x, p) plane for 9 E [0, 2tt). The first 
order equations are 



dx 
dx 
dx 
dp 

dx 



1 - 02 _ r{x^ + p2) 

— 2r0x 
— 2r0p 



(45) 
(46) 
(47) 



There are solutions connecting the minima (±1,0,0). 
One of them is given by x = P = and 4'{x) = tanh(a;). 
Its energy is 4/3, and it may give rise to a domain wall in 
(3, 1) space-time dimensions. There are other solutions 
connecting the minima (±1,0,0), with the same energy 
of the above solution. They are given by 



p{x) 



tanh(2ra;) 



± 



± 



cosh(2rx) 

sin6' 
cosh(2rx) 



(48) 
(49) 

(50) 



They form an elliptical surface, and in Fig. [2] we depict 
these solutions in the case r = 1/6 and 6 = 7r/6. 
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FIG. 2: The BPS solution described by Eqs. 
@, for r = 1/6 and 61 = 7r/6. 



), (toh, and 



We notice that this three-field solution may be used 
to describe three-wave parametric solitons in quadratic 
nonlinear media |2^]. Moreover, the fact that x and p are 
non- vanishing for a; — > indicates that the wall generated 
by the (f) field may entrap both x and p in its interior -- 
see Refs. [^2[ ^ ^ for further details. 

There are other solutions, connecting the minima 
(±1,0,0) to (0, -y/l/rcos^, ^/T/rsm9) through parabolic 
surfaces. Some of them are, for r = 1/4, and for 

92^(l/4)(l-0) 



X 



4>{x) 



1 



-(1 + tanh(a;/2)) 



X{x) = ± cos 9yj2{\ - tanh(a;/2)) 
p{x) = ± sin 6*^2(1 - tanh(j:/2)) 

and also, for x^ + ^ (1/4) (1 + 0) 



(x) = --(l-tanh(a;/2)) 



X{x) = ± cos 9yj2{\ + tanh(a;/2)) 
p{x) = ± sin 9yj2{l + tanh(x/2)) 



(51) 

(52) 
(53) 

(54) 

(55) 
(56) 



These solutions have the same energy, given by the value 
2/3. This is one-half the value of the energy in the former 
sector, defined by the minima (±1,0,0). 

We notice that the potential generated by T4^(0, x, p) 
as in Eq. (H) is such that 



V{Q,X,p) = \r^i\ 



X 



(57) 



It shows that the domain wall generated by the (f) field 
may entrap the other two fields in its interior, giving rise 
to a model where one can have SSB of the continuum, 
global U{1) symmetry. In Eqs. (|45|), (Ei), and (li^), if we 



get to the equations d(j)/dx — 1 



and da/dx 



—2r(f>a, which are similar to the first order equations that 
appear in the model represented by Eq. (22), thus we can 
follow Ref. j25j to find all the BPS states of the model. 

The continuum, global U{1) symmetry that appear for 
s = can be further explored. For instance, we notice 
that in (3, 1) space-time dimensions the host wall (t>{x) 
depends only on x, and may entrap the other two fields 
in two spatial dimensions. Now, if we let the continuum 
symmetry that appear in the (%, p) plane to become local, 
and if we start with all the needed physical ingredients, 
the model gives rise to a host domain wall that entraps 
a planar system, driving us for instance to the scenario 
where the so-called Callan-Harvey effect |^ may spring 
very naturally. We recall that the Callan-Harvey effect 
concerns anomaly cancelation in (2, 1) space-time dimen- 
sions due to the presence of a domain wall mass to a 
fermionic field, and we quote Ref. |Q for an Abelian 
version of this effect. This issue is presently under inves- 
tigation, and we hope to report on it in the near future. 

Another line of investigation concerns the high tem- 
perature effects, which may drive the system to phase 
transitions. This line may follow the two last work in 
Ref. , which deal with the thermal effects in the case 
of two real scalar fields. Furthermore, since the high tem- 
perature effects are obtained a la Matsubara, which is a 
compactification process, we can also consider more gen- 
eral possibilities in which one investigates the effec- 
tive potential in the case two or more dimensions of the 
original space-time are compactified. Similar ideas have 
been introduced in |29 for bosons and fermions, examin- 
ing the temperature inversion symmetry. This direction 
brings other issues since together with the usual thermal 
effects, we are also investigating effects of compactifica- 
tion of one or more spatial dimensions. 

The presence of topological sectors depends on the 
manifold of vacuum states of the system, and this indi- 
cates that our investigations are model-dependent. Thus 
we can build other three-field models, exploring the pos- 
sibility of finding spatial junctions of defects, which may 
allow that we fill space with solid structures, generalizing 
the idea of tiling the plane with regular poligons. In the 
plane there are an infinity of regular polygons, although 
only three of them are capable of tiling the plane. In 
space, however, there are only five regular polyhedra - 
the five Platonic solids - and the issue is then which are 
capable of filling space, and how we can generate them 
within the context of three-field models. These and other 
related issues are presently under consideration, and we 
hope to report on them in the near future. 

We would like to thank A.R. Gomes, J.R.S. Nasci- 
mento and R.F. Ribeiro for discussions, and CAPES, 
CNPq, PROCAD and PRONEX for partial support. 
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